Geometric Correlations and Breakdown of Mesoscopic Universality in Spin Transport 
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We construct a unified semiclassical theory of charge and spin transport in chaotic ballistic and 
disordered diffusive mesoscopic systems with spin-orbit interaction. Neglecting dynamic effects of 
spin-orbit interaction, we reproduce the random matrix theory results that the spin conductance 
fluctuates universally around zero average. Incorporating these effects in the theory, we show that 
geometric correlations generate finite average spin conductances, but that they do not affect the 
charge conductance to leading order. The theory, which is confirmed by numerical transport calcu- 
lations, allows us to investigate the entire range from the weak to the previously unexplored strong 
spin-orbit regime, where the spin rotation time is shorter than the momentum relaxation time. 

PACS numbers: 73.23.-b, 85.75.-d, 72.25. Dc 



At low temperatures, linear electric transport proper- 
ties of complex mesoscopic systems are statistically deter- 
mined by the presence of few symmetries only, most no- 
tably time- reversal and spin rotational symmetry 0, Q • 
This character of universality is believed to be indepen- 
dent of the source of scattering in the system, and to exist 
in both ballistic chaotic quantum dots or diffusive disor- 
dered conductors Q. Universality in electric transport 
holds not only for global properties such as the conduc- 
tance, but also for correlators of transmission amplitudes 
between individual channels. Thus, it is natural to ex- 
pect that all transport properties that depend solely on 
the scattering matrix are universal as well. This conjec- 
ture has been theoretically verified for all charge trans- 
port properties, under the sole assumption that scatter- 
ing generates complete ergodicity. Inspired by Ref. 
several recent theoretical works [B|~0] have further sug- 
gested that spin transport in mesoscopic systems with 
spin-orbit interaction (SOI) also displays universal ran- 
dom matrix theory (RMT) behavior. The agreement be- 
tween numerics for a disordered lattice H and the RMT 
prediction Q| for the mesoscopic fluctuations of the spin 
Hall conductance indeed seems to imply that RMT uni- 
versality also exists in magnetoelectric transport. 

In this work, we go beyond the conventional semiclas- 
sical theory of transport and show that even when all 
requirements for universality are met and the fluctua- 
tions of the spin and charge conductance as well as av- 
erage charge conductance remain universal, the average 
spin conductance is finite in disagreement with the RMT 
prediction. This effect originates from the spin-orbit cou- 
pling through which the electron spin perturbs the elec- 
tron dynamics in such a fashion that, certain dynamical 
correlations survive despite the self-averaging nature of 
ergodic dynamics. These correlations depend on the ge- 
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Figure 1: (Color online) Spin-dependent transmission coeffi- 
cients 7%i, Eq. Q, for (a) weak and (b) extended range of SO 
coupling showing the crossover from cubic (green line) to lin- 
ear (blue line) behavior for the two-terminal chaotic quantum 
dot shown in the inset of panel (a), (c) Spin current polar- 
ization angle 9 = arctan(Tfi /Tli ), for the system shown in 
the inset, where the right lead encloses an angle of <f> with the 
z-axis in the linear (cubes) and the cubic (circles) regimes. 



ometry of the system, namely the relative positions of 
the leads connecting the system to external electronic 
reservoirs as well as the form of the SOI. As an exam- 
ple, we consider a two-dimensional quantum dot with 
Rashba SOI and find that the average two-terminal 
spin-conductance G' 1 is proportional to (z x R) M . Here 
the vector R connects the two terminals, z is the unit 
vector perpendicular to the dot and fx is the spin com- 
ponent. This is illustrated in Fig. Ufa) for the corre- 
sponding spin transmission. The polarization of the av- 
erage spin current is thus determined by the direction 
of the average electronic flow. In bulk diffusive systems, 
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when the mean free path is shorter than the spin rotation 
length, this effect reduces to the extraction of the current- 
induced spin accumulation (CIS A) and the spin Hall ef- 
fect fiol- 12 1 in finite systems. We stress however that the 
consequences of these geometric correlations have been 
considered in neither charge nor spin-transport in quan- 
tum dots. Moreover, our calculations extend the exist- 
ing theory for CISA and SHE in finite diffusive systems 
to the strong SOI regime (i.e. mean free path is longer 
than the spin rotation length). It is of practical impor- 
tance to point out that the process that leads to finite 
spin conductance is robust against temperature smear- 
ing or dephasing. From the point of view of mesoscopic 
spintronics, this opens up possibilities towards an elec- 
trically controlled generation and detection of pure spin 
currents, since the uncontrolled mesoscopic fluctuations 
are suppressed by simply raising the temperature. 

We consider a mesoscopic quantum dot with no par- 
ticular spatial symmetry as sketched in the insets of 
Fig. [T] We treat impurity and boundary scattering on 
equal footing and consider diffusive as well as ballistic 
chaotic charge dynamics. The dot is connected to two or 
more external leads. For simplicity, we assume idealized 
reflectionless leads in which the SOI vanishes. The real- 
istic case of finite SOI in the leads can then be obtained 
by combining the scattering matrices of the realistic leads 
with that of the quantum dot. This choice allows us to 
uniquely define transport spin currents through a cross- 
section of the leads without the ambiguities that plagued 
bulk calculations [3] . The leads are maintained at differ- 
ent electrochemical potentials eVi, but have no spin ac- 
cumulation. The scattering approach to transport gives 
the spin and charge currents in lead i as [TjJ 



e 2 

3 



\Vi-Vi), 



(1) 



with the generalized, spin-dependent transmission coeffi- 
cients 7?? obtained by summing over all transport chan- 
nels in leads i and j ' [J, 0| , 
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H,u = 0,x,y,z. (2) 



Here, are Pauli matrices (ero is the identity matrix) 
and the trace is taken over the spin degree of freedom. 
The transmission amplitudes in Eq. ([2]) can be expressed 
in terms of the Green's function |15j . Next, we obtain the 
full Green's function G R (r, r') by either (i) the conven- 
tional Born approximation for impurity scattering inside 
the conductor or (ii) by a multiple reflection expansion 
for boundary scattering case (ii), G R (r, r') is 

expressed as an iterative solution of 

G R (r,r') = G R (r,r')-2 [dadG R (r, a)G R {a, r'), (3) 



we evaluate the surface integrals in Eq. ^ asymptoti- 
cally as kpL oo, where kp is the Fermi wavenumber 
and L is the linear size of the conductor (lTjl . We obtain 
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where the sums run over all trajectories 7 starting at yo 
on a cross-section of the injection lead and ending at y 
on the exit lead. The classical action of 7 is 5 7 in units 
of h and its stability is given by A~ j which includes a 
prefactor (27rifi,) -1 / 2 as well as Maslov indices. For the 
spin dependent part, we specialize to the Rashba SOI 
Hji = (hk a /m)(p x a y —p y <J x ), where fc" 1 is the spin pre- 
cession length [9|. We then obtain 



v i = n v ^ = n + sUi >-y + 1 5u t 

i=l i=l 

k(^ { &in(fcfj r. 



hw\ 

7 1 



SUi 7 = — — 



4fc F 

k a j sin(fc Q Fj 
2fc 



- 1 ] rj ■ Vi 
1 I 



cos 6j 



cr, sm ( 



2fcplr.il cos 6, 



(1 - cos(fc Q |ri|)) . 



(6) 
(7) 



Here £ = for a disordered system with weak, short- 
ranged impurities and £ = 1 for a ballistic quantum dot 
with hard-wall confinement or a disordered system with 
strong, extended impurities. In both cases 7 consists of 
segments r. t = (xi,y i: 0) with % = 1,2, ..JV 7 , f< = r,/|rj|, 
n, is the (inner) unit normal vector and Oi is the an- 
gle of incidence at the ith reflection point, rj = z x cr 
and Ui.j = exp[— ik a r) ■ is the Rashba spin rota- 

tion matrix along that segment. We note that there are 
also corrections to A 1 which we have already ignored 
here, because they do not contribute to the spin conduc- 
tance (however they generate diffractive corrections to 
the charge conductance). The Eqs. (j4][7|) completely de- 
scribe spin and charge dynamics of coherent conductors. 

The conventional semiclassical theory is obtained via 
the approximation V 1 ks ]~Ii=i ^vy> which leads to the 
universal RMT predictions for charge transport 3, HI- 



We now show that this approximation also leads to RMT 
results for spin transport for fj, =/= 0. We first start from 
the diagonal approximation, where 7 = 7', and obtain 
tr^cT^V^'] = 0, showing that the diagonal contribu- 
tion to the spin current vanishes. The next-order con- 
tributions within the conventional semiclassical theory 
of transport are the loop corrections, in which a self- 
crossing trajectory 7, is paired with a path 7' avoid- 
joing around the loop in the the 



ing the crossing and 
opposite direction [2C 
time-reversed of 7, and the loop contributions are pro- 



2l|. Along the loop, 7' is the 



where dG R (r, a) — h a ■ VG R (r, x)| x=Q , with n Q the (in- portional to (tr[U 7l a^U 7l ]) 
ner) unit normal vector at the boundary point a. Finally, tation along the loop only. For large SOI, C/ 7i is random, 



where E/ 7 , gives the spin ro- 
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thus averaging produces vanishing weak localization cor- 
rection to the spin conductance. For weaker SOI, we 
expand all spin rotation angles to second order in k a L to 
obtain (tr[J7 7i cr M J7 7i ]) ~ 2i6 MZ ( sin (fc 2 &A 7i )) . The area 
difference 5A~ ll is given approximately by twice the di- 
rected area of the weak localization loop. For a chaotic 
system, the areas are symmetrically distributed around 
zero, thus the average vanishes. We note that extending 
the semiclassical approach of Ref. [22| to the calculation 
of the variance of the spin conductance, one straight- 
forwardly reproduces the leading-order RMT results of 
Ref. [H . Details of this calculation will be presented else- 
where 1 17). We conclude that conventional semiclassical 
theory, which neglects effects of spin on the charge dy- 
namics, only reproduces RMT predictions. 

We next include the effects of SOI on the electronic dy- 
namics and consider a two-dimensional conductor which 
can be either a ballistic quantum dot with hard-wall con- 
finement, or a disordered system with short-ranged im- 
purities. To do this, we go back to Eqs. ((IE]) and include 
the corrections to the amplitude A and the spin matrix 
U to order 0{k a /k-p) and 0(1/ /cf l r i I ) ■ We now assume 
that different trajectory segments are uncorrelated and 
define Uz j7 = nS+i ^vr t° obtain 
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We see that spin currents have contributions from ev- 
ery trajectory segment, which are further rotated by 
the fluctuating spin-orbit fields of the subsequent reflec- 
tions. We distinguish three different regimes that de- 
pend on the balance between linear system size L, the 
mean distance between (boundary or impurity) scatter- 
ings £ = (|r/|), and SOI length A:" 1 : (i) the spin-ballistic 
small SOI limit k a L, k a £ -c 1, (ii) the spin-diffusive limit 
k a £ < 1 <C k a L, (iii) the spin-chaotic strong SOI limit 
k a L,k a £ ^> 1. In regimes (i) and (iii), the orbital dy- 
namics can be chaotic ballistic or diffusive depending on 
the ratio between £ and L. We will be focusing on long 
ergodic or diffusive trajectories 7 for which we ignore the 
averages (sin0j) 7 and (n,) 7 for all three regimes, save for 
the case of a quantum dot in regime (i) (see below). 

In the small SOI regime (i), we expand the rhs of 
Eq. ([8]) to leading order in k a £ setting = 1 and 
1 - sin(fc Q |r i |)/fc Q |r i | ~ (k a \n\) 2 /6 in Eqs. (JSC]). We get 
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We now perform the average (. . .) 7 over the set of tra- 
jectories 7. Although individual r.; are pseudorandom 
in length and direction, being generated by the cav- 
ity's chaotic dynamics, they satisfy J2i r l ~ where 
Rjj is the 7-independent vector connecting the injec- 
tion and exit terminal. We thus obtain (tr[V^ <t m ]) = 



C[kl£(l + 2£)/(3fc P )](z x R„) M , where C is a number of 
order one that depends on geometric details of the cav- 
ity. This factor multiplies the independently averaged 
orbital terms in Eq. ^ for 7 = 7', which we compute as 
in, e.g. Ref. [2l|. We estimate £ = (|r ; |) ~ irA/L for a 
chaotic dot of area A and perimeter L, and £ = vpr for a 
diffusive system with momentum relaxation time r. We 
finally obtain 




i«L. 



(10) 



with the number N{ = Int(fcpWi/7r) of channels in lead i, 
TVt = Ni and W = minW^ the width of the narrowest 
lead. In the ballistic limit, this formula has a correction 



term 



3k F N'i, 



J2i Ni(z x R) M , where R; is the average 
momentum direction of electrons entering through lead 
I, originating from nonzero (fij) 7 [17|. We see that the 
average spin-dependent transmission, and thus the aver- 
age spin currents, are determined by the relative position 
of the injection and exit lead and are proportional to the 
classical conductance from j to i. 

In the spin-diffusive case (ii), L 3> k~ 1 3> £, the spins 
precess around randomly oriented SOI fields, thus relax- 
ing via the Dyakonov-Perel mechanism. In particular, we 
can no longer set 1X/ j7 = 1 in Eq. ([5]). Instead, we assume 
that 7 is a stochastic sequence of segments with random 
orientations tpi, which determine the spin rotation E/j j7 . 
The sequence of rotations is computed by averaging over 
tpi . For a general Pauli spin matrix s • cr one has 

J ^U^s-aU}^ = cos 2 (fc Q |r ? ;|/2)s.<T (11) 

^|r l | 2 /2)sin 2 (fc Q |r ? ;|/2) J7 (s-cr)r 7 . 

This average is different for in-plane and out-of-plane po- 
larization, which is the origin of the anisotropy of the 
Dyakonov-Perel spin-relaxation time. In our case, the 
generated spin is in-plane and the second term in Eq. (|11[) 
vanishes |24|. We have 



l + 2£ 




V ■ r i 



(12) 



where we used k a £ <C 1 <C k a L, approximated |r^| « I, 
Vi and introduced the duration r; of the first I segments 
of 7. For each possible choice of I, the spin rotation 
thus separates into a spin independent piece for segments 
1, • • ■ ,1 — 1, a spin generation piece on segment I, and a 
spin relaxation piece on segments I + 1,- •• , 7V 7 . Fix- 
ing the endpoint r/ of segment I and summing over all 
possible orbits we obtain that the spin conductance is 
proportional to a product of (i) a diffusive probability 
P(xi,Xj) to go from the injection lead to r;, (ii) a spin 
generation factor (1 + 2^)fc^ £77 • (x; — x;<)/12fcp multiply- 
ing the probability of ballistic propagation from x; to , 
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(iii) a diffusive probability to propagate from point x;/ to 
the exit lead times the probability that the spin survives 
this diffusion. Thus we have 

k 3 £ f 

(7^°) oc £ 3/^-^7 / dx i dx.jdxidx il P(x i ,Xj){Ki--x 1 >)' / 

p -\xi-x t /\/i 

x (1 + 2^— f P(x 4 ,x i ,)e- fc °l x "- x 'l. (13) 

Since the length scale characterizing P(x/,Xj) is L, we 
evaluate the integrals above asymptotically in the limit 
k a £ «C 1 <C k a L. After some algebra we finally obtain 



cx sgn(fc Q )(l 



^(^x^), (14) 



up to a factor of order unity depending on details of how 
the leads (with width W) are attached to the cavity. Not- 
ing that for our geometry is in the direction of the 
current flow and its magnitude is L, we obtain that the 
spin conductivity is a s oc ek\£ 2 in agreement with the 
spin diffusion equation calculations Icj, 12 1 . 

Spin chaos regime (iii): Similar to regime (ii). we 
average over uncorrelated direction angles 9i but do 
not Taylor-expand sin(fc a ,|r;|)/fc a |ri| — 1. We instead 
take the average over the segment lengths |r|,; as 
n^J ;+1 (cos 2 (fc Q |ri|/2)) « 2 N ~<~ 1 in a chaotic/stochastic 
system with k a L 1. Eq. (fT2"|) is then replaced by 



1 + 2$ 



k a ( sin(fc a |r;| 
2fc F I fcJrjl 



- 1 p7 • r i 



Averaging over 7 we see that the dominant contribu- 
tion is the last term. We thus approximate the sum by 
its last term, and take fc Q |rjv_,| ~ k a L 3> 1 to obtain 
{Yy K y t } 7 = l + {C'k a /2k F )r] ■ R r Here C is (1 + 2£) 
times a constant of order unity that depends on the de- 
tails of the scattering near the lead. We finally obtain 
the transmission coefficient 



k F W£/L £<L. [ 0) 



Equations (fTU]) . (jT5J) and (fT5|) are our main results. 
They show how a finite spin conductance emerges from 
classical geometric correlations depending on the posi- 
tions of the leads. These equations can be straight- 
forwardly extended to Dresselhaus SOI by substituting 
z x Q -> (Q x , -Q y , 0) for Q = R„ [Eqs. PJ) and fT3])] 
or Q = Rj [Eq. ijlgjl]. 

To check these predictions we performed numerical re- 
cursive Green's function quantum transport calculations 
for a tight-binding Hamiltonian [23[ with Rashba SOI 
and evaluated the spin-resolved transmission probabil- 
ities between two leads as defined in Eq. ([2]) for both 
the chaotic and diffusive cases. We computed the trans- 
mission for chaotic cavities, shown as insets in Fig. [TJ 



averaged over 2000 different configurations of the Fermi 
energy and the position and orientation of the central 
antidot. Panel a) shows for the small a = ak a regime (i) 
that the numerically obtained (dots) for the cavity 
in the inset agrees very well with the predicted cubic be- 
havior, Eq. (fTU)) , (solid line) for C = l. In panel b) Tf° is 
depicted for the same chaotic cavity (black circles) and 
for a square cavity with Anderson disorder (violet trian- 
gles) for the entire range from weak to strong SOI (regime 
(i) to (iii)) demonstrating the crossover from cubic to lin- 
ear behaviour according to Eqs. (fTUj) and p3|) . In panel 
c) we further numerically confirm the predicted direction 
of the in-plane spin polarization 6 = arctan(T| 1/^2 1) f° r 
regime (i) (dashed line, Eq. (fT0|) ) and regime (iii) (solid 
line, Eq. (fT5|) ) by rotating the right lead around the semi- 
circle billiard shown in the inset. 

In conclusion, we have presented a semiclassical calcu- 
lation of spin transport in mesoscopic conductors which 
incorporates next-to-leading order corrections to the 
semiclassical Green's function. We showed that in con- 
trast to RMT predictions, the average spin conductance 
does not vanish, even if all the conventional conditions 
for universality are met. Our method moreover allowed 
us to investigate the strong SOI regime for finite diffusive 
systems for the first time, Eq. (fTSj) . 
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